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Abstract 

We study the systematic inclusion of QED corrections in the evolution of parton 
distributions. 0{a) corrections modify the evolution equation for parton distribu- 
tions. They introduce additional parton distributions, like the photon distribution 
in the nucleon, and lead to additional mixing effects. We discuss the modifications 
for a realistic model of A^/,up up-type flavours and Nf ^down down-type flavours. We 
have implemented these corrections into a numerical program and we quantify the 
size of these effects in a toy model. The corrections reach the order of 1%. 



1 Introduction 



Electroweak radiative corrections to QCD observables can be sizeable and have to 
be taken into account for precision measurements. If the observable involves hadrons 
in the initial state, this does not only imply the inclusion of corrections to the hard 
scattering process, but also to the parton distribution functions. Although these parton 
distribution functions parameterize the non-perturbative information, they depend on a 
factorization scale. This dependence on that scale can be calculated within perturbation 
theory. Usually parameterizations of parton distributions functions are given at a low 
input scale and are then evolved to a suitable scale for the considered scattering process. 

The inclusion of QED corrections modifies the evolution equations. Although this 
issue has already been partially addressed in the literature [1,2], a complete and detailed 
discussion is still missing. With this letter we hope to fill this gap. The inclusion of 
QED effects leads to additional partons in the nucleons: Quarks have an electric charge 
and may radiate off photons. Photons in turn may split into a pair of charged leptons. 
Therefore, the set of partons inside the nucleon is enlarged by the photon and the charged 
leptons and anti-leptons. 

In pure QCD the evolution equations are a highly coupled set of equations, if one 
chooses as a basis the gluon, quark and anti-quark distributions. From a practical point 
of view it is advantageous to choose a different set as a basis for the parton distributions, 
such that most distributions evolve multiplicatively and mixing occTirs only in a minimal 
sector. In pure QCD such a set is provided by the gluon distribution, the quark singlet 
distribution and the remaining quark non-singlet distributions. Mixing occurs here only 
between the quark singlet and the gluon distribution. The inclusion of QED effects leads 
to additional mixing and we discuss in detail the modifications for a realistic model of 
A^/,up up-type flavours and iV/down down-type flavours. We show that the inclusion of 
QED effects leads to mixing between the quark singlet, the gluon, the distribution A[//3, 
being the difference of all up-type flavours minus all down-type flavours, the photon, and 
the distribution E;, being the sum of all charged leptons. It is worth to summarize the 
sources of all modifications due to the inclusion of QED effects: 

• The electromagnetic interaction enlarges the set of partons inside the proton to 
include the photon together with the charged leptons and anti-leptons. 

• Up- and down-type quarks have different electrical charges. This is in contrast to 
the pure QCD case where all fermions carry equal colour charges. The difference in 
the electrical charge leads to additional mixing compared to a hypothetical model, 
where all quarks would have equal electrical charges. 

• The number of up-type quarks iV/,up does not necessarily equal the number of down- 
type quarks A^/,down- Again, a non-equal number of up- and down-type quarks will 
lead to additional terms in the evolution kernel as compared to a model with an 
equal number of up- and down-type quarks. 

We have implemented the QED corrections into two numerical programs for the evolution 
of parton distributions. One is based on program described in Ref . [3] . We use the A'"-space 
method with a parabolic contour [4] . 
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We work, as it is usually the case for QCD processes, with the MS factorization 
scheme. The inclusion of QED corrections to an observable with hadrons in the initial 
state requires the QED corrections to the hard scattering process as well as to the evolution 
of the parton densities. We would like to point out, that for correct results both corrections 
have to be calculated in the same scheme. 

This paper is organized as follows: In the next section, we shortly review the theoretical 
background and we present the evolution equations with the inclusion of QED effects. In 
Section |21 we quantify the size of these effects in a toy model. Finally, Section |^ contains 
a summary and the conclusions. 

2 Parton evolution including QCD and QED corrections 

The evolution equation for the parton distributions /(x, Q^) of the proton reads 

Q'^-^^l0^ = PM')®fM% (1) 

where x stands for the momentum fraction of the nucleon carried by the parton, P{x, Q"^) 
is the Altarelli-Parisi evolution kernel, and ® denotes the convolution 

1 1 

(A O B){x) = J dy J dz 5{x - yz)A{y)B{z). (2) 



In general, /(x, Q"^) is a vector consisting of M different parton distributions and P{x, Q^) 
is a M X M matrix. In order to study QED effects in the evolution equation, we have 
to enlarge the usual set of partons consisting of light quarks {u,d, s,c,b), anti-quarks 
{u,d,s,c,b), and the gluon (g) by charged leptons (e,/i,r), charged anti-leptons {e,Ji,f), 
and the photon (7). 

We define the valence distributions by 

/. = / - /, (3) 

where / is an arbitrary fermion. In pure QCD with 5 flavours it is convenient to take the 
valence distributions together with the flavour- singlet distribution 

s= E (^ + ^1 (4) 

q=u,d,s,c,b 

and the following non-singlet combinations 

A2 = u + u — d — d, 

A3 = u + u + d + d-2{s + s), 

A4 = u + u + d + d + s + s — 3{c + c), 

= u + u + d + d + s + s + c + c-4:(b + bj (5) 

as a basic set of parton distributions. The inclusion of charged leptons and the photon 
enlarges this set and leads to the basis 

Basis 1 = {'^t;? du^ s^^ c^, 6^, A2,A3,A4,A5,E A'2,Af5,E',7} 
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with the leptonic quantities 

A2 = e + e — fi — Ji, 

A?5 = e + e + /i + /i-2(r + f), 

T}= y: 0+0- (6) 

This basis has the advantage that in absence of QED effects all parton distributions 
except the one for the quark singlet combination E and the gluon distribution g evolve 
multiplicatively. Mixing only occurs between E and g. However, in this basis the inclusion 
of the QED effects leads to mixing between eight parton distributions, namely A2, A3, 
A4, A5, E, g, El and 7. 

To avoid the evaluation of an 8 x 8 matrix in the evolution kernel, we use a different 
basis. 

Basis 2 = {uy, d^, s^, c^,, 6^, A^^, A„c, A^^, Aud, g, e^, /i^,, r^, A2, A3, E', 7}, 

which leads to a mixing between only five distributions. The new quark combinations in 
Basis 2 are defined by 

A(7D = u + u + c + c — d — d — s — s — b — b, 

= u + u — c — c, 
Ads = d + d — s — s, 

Asb = s + s-b-b. (7) 

AuD is the difference between all up- type and all down- type quarks, while Auc, Ads and 
Asb parameterize the differences between two up-type quarks or between two down-type 
quarks. This basis takes the different electric charges of up- and down-type quarks into 
account and leads to mixing among only Aud, E, g, 7 and E/.^ 
Using Mellin moments 

1 



r = /dxx^-V(x), 



the evolution equation (^Q) factorizes into the form 

Q^^^^ = p^(g^)-r(g^). (9) 

In this work, we consider the evolution kernel to next-to-leading order (NLO) in = 
as/ 4:7!- and to leading order in a = a/ Arc, 

P'^m = a.(g')Po^QCD + + «W')^(fQED + 0{al) + 0(aa.) + Oia'). 

(10) 



^It should be noted that Basis 2 has the sHght disadvantage that in the absence of QED effects and 
for an unequal number of active up- and down- type flavours it leads to mixing between Aud, E and g. 
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Mixed terms of order 0{aas) are not included. We separate the kernel into a QCD piece 
and a QED piece, 



^qed(Q') = a(Q')^o'iED + 0{a'). (11) 

The QCD evolution kernels -Pq^cd ^^"^ A^cd ^^^^ known [5-13]. Explicit expressions 
for the anomalous dimensions contributing to -Pqqcd ^"^^ A^cd can be found in Ref . [9] . 
The QED evolution kernel PqIqed is easily obtained from PqIqcb by adjusting the colour 
factors, 

pN _ f^2 pN pN _ r\2 ]\rf pN 

-Po)qED,7/ = <5/-Po,7/' -Po^QED,77= T, [ Q) Po,^-f i (12) 

where Q/ is the relative electric charge of the fermion / and N^. denotes the multiplicity 
due to colour degrees of freedom, i.e. N[ = 1 for leptons and A^/ = 3 for quarks. The 
QED evolution kernels without prefactors read 



pN 



45i(iV) - 3- 



N{N+1) 



N _ o N^+N+2 



0,/7 ^Af(Af+l)(iV+2) ■ 



pN _ 9_/V^+JV+2_ pN _ _4 /iqx 

-^0,7/ — ^(Ar-l)Ar(Ar+l) ' "^0,77— 3 

with Si{N) =Y.n=l^/n. 

The evolution equations of the couplings and a are given by 

Q^^ = /5K), Q^^ = /5QED(a), (14) 

where the beta functions for the strong and electromagnetic couplings are expanded to 
the appropriate order, 

/3(a,) = -/5oa2-/5ia^ + 0(at), 
/?QED(a) = -/^o.qedo^ + C(a=^) (15) 

with the lowest-order QED beta function 

/?0,QED = -^ E ^IQ)- (16) 

f=e,^,T,u,d,s,c,b 

In pure QCD there are several prescriptions, how to truncate the perturbative expansion 
for the evolution kernel and the beta function. Although they all agree formally to a 
fixed order in perturbation theory, they differ numerically due to a different treatment of 
higher order terms. In this paper we use the direct solution of Eq. together with an 
approximate solution for a^. In Ref. [3] this method is labelled "x-space and truncated 
in 1/L solution". 

Using the operator 

E^(Q^Q^) =rQ.exp j dtP^(e*), (17) 

InQg 
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the evolution equation ^ is immediately solved by 

nQ') = E''iQ\Ql)nQl). (18) 

The symbol Tq2 in ()17p means ordering in Q^. The evolved parton distributions in x-space 
are then obtained by an inverse Mellin transformation, 

p+ioo 

f(x,Q') = ^ [ diVx-^r(g2), (19) 

p—\oo 

where p has to be chosen such that all singularities lie on the left-hand side of the inte- 
gration contour in the complex plane. 

The valence distributions as well as the distributions A^^, A^i^, ^sb-, ^2 and A3 
evolve multiplicatively. For the and distributions we have 

Q'^d^^iQ") = KcD,n,,=-i(Q') + a{Q')QlPo'!ff) ■ d^HQ'), (20) 

where -PQCD,ns,r;=-i(Q^) is the NLO QCD evolution kernel for a non-singlet valence distri- 
bution, 

PQCD,ns,ri=~liQ ) = (^s{Q )-Po,QCD,ns + '^s (Q ) -Pl,QCD,ns,»y=-l' (21) 

and rj = —1 selects the odd moments corresponding to valence-type quark distributions 
in -Pi^QCD.ns- QED evolution kernel, we explicitly factored out the coupling a{Q^) 

and the electric charge in (|20|) . Similarly, the evolution equation for A^^, Auc and Asb 
read 

Q'^^AUQ') = {pScu,ns,,=im+amQiPo'!ff) ■ ^im, 

Q'^^uciQ') = (^QCD,„s,,=i(Q') + aiQ')QlPo'!ff) ■ Al{Q'), 
Q'^2^sb{Q') = {PScB,ns,,=i{Q') + amQlP^^jf) ■ Af,(g^). (22) 

Here t] = 1 selects the even moments corresponding to non-singlet non- valence-type quark 
distributions. The remaining non-singlet evolution equations can be obtained in a similar 
way. For the evaluation of (fTH|) we closely follow Refs. [3,4]. In particular, we use the QCD 
evolution operator defined in Eq. (8.26) of Ref. [4] for the quark non-singlet distributions. 
The QED part of the evolution operator factorizes in this case and we have 

pN 
O.QED 

In the rest of this paper we set the valence distributions Sy, Cy, by, Cy, fiy and r„ as well 
as the lepton combinations A2,/ and A3 ; to zero. 
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As already mentioned, the inclusion of QED effects lead to a mixing among the parton 
distributions Ajj^), E, g, 7 and S'. This is in contrast to the pure QCD evolution, where 
mixing occurs only between the singlet combination E and the gluon distribution g. It is 
convenient to separate the evolution kernel into a QCD and a QED piece as done in (fTT|) . 
Let f^{Q^) be the vector 



N 



(24) 



The QCD part has the form 



/ pN 



pN 



pN 

^QCD,qq 



pN 
QCD,ns,?y= 



pN 
^QCD,qq 
pN 

^QCD,qq 






^NfPQCD,qg 0^ 


0/ 



(25) 



where 



/,UP 



/,down 



/jdown 



(26) 



and PQCD,qq, PQCB,qg, ^QCD.gg aud PqcD.gg are the components for the usual evolution 
matrix in the basis {^,g} without QED corrections. Note that the difference Pqqd qq ~ 
^QCD,ns,,7=i is of Order al 

For the QED part it is convenient to introduce the sum and the difference of the 
squares of the electrical charges of up- and down-type quarks. 



\{Ql 



Q 



(27) 



Using these definitions, the QED evolution kernel Pq^qeb yields 



0,QED 



fc^P^ 



ojf 

pN 



pN 
-^0,7/ 



V 



.P.^jf 2N,Nf 













Nf 
+ C 



+ C_ 



pN 
pN 



3 /o pi 
— 4P0,QED-r(), 



.77 



2N,P, 



0\ 



pN 

pN I 



:28) 



where iV/ is the number of charged leptons. Note that for c_ = and bjq^ = 0, e.g. equal 
charges and equal number of up- and down- type quarks, decouples and evolves 

multiplicatively. It is worth to point out that c_ parameterizes the difference in the 
electrical charges of up- and down-type quarks, and that 6^^. parameterizes the difference 
between the number of up-type quarks and the number of down-type quarks. Both effects 
introduce additional mixing terms, as can be seen from Eq. (j28j) . 

For the numerical evaluation, we have written two independent programs which are 
in numerical agreement. 
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3 Numerical results 



For the numerical discussion, we adopt the toy model of Ref. [14] which consists of 
four active flavours with Aqcd = 250 MeV. At the input scale Qq = 2 GeV, the initial 
parton distributions are as follows 



X 



c S = Asx'^ ^l - xY, xg= AgX-'^-^l - xf, 

xc=0, xc=0. (29) 



The sea S is taken to be SU(3)flavour-symmetric and carries 15% of the nucleon momentum 
at the input scale. These assumptions, together with the usual flavour and momentum 
sum rules, fixes the constants Au, Ad, As and Ag. The sea is related to the quark singlet 
distribution by 

S = - Uy - dy, s = (q + q) ■ (30) 

q=u,d,c,s 

Furthermore, SU(3)flavour-symmetry of the sea S implies 

u = d = s (31) 
at the input scale Qq. For the running of the strong coupling we use 

To study QED effects, we introduce some further assumptions. We require that 7 and 
vanish at the input scale Qq, 

7 = 0, = ^ (/ + r) = 0. (33) 

l=e,fi,T 

For the running of the electromagnetic coupling we use 



a{Q') = 775n (34) 

l + a(m2)/3o,QEDln(5) 

with 

a(m2) = 1/133.4, = 1.777 GeV. (35) 

Since this model has only four active flavours, all bottom quarks have to be omitted in 
((Zj) and therefore by and Agt are ignored in the calculation. The lowest-order QED beta 
function in this case takes the value 

/?0,QED = (36) 

To verify the numerical precision of our evolution code, we have checked, that the nu- 
merical values in the absence of QED corrections, but with a 5 x 5 mixing matrix as in 
Eq. (j25j) . agree with the benchmark values quoted in Ref. [14] as well as with the values 
obtained from the numerical program in Ref. [3] . 
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Figure 1: Difference in per cent due to QED effects of the w-quark, d-quark and gluon 
parton distributions as a function of x for Q = 100 GeV (l.fi.s.) and as a function Q for 
X = 0.1 (r.li.s.). 
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Table 1: Differences in per cent for tlie evolution from Qq = 2 GeV to Q = 10 GeV due 
to the inclusion of QED effects. 
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Table 2: Differences in per cent for the evolution from Qo = 2 GeV to Q = 100 GeV due 
to the inclusion of QED effects. 
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Table 3: Differences in per cent for the evolution from Qq = 2 GeV to Q = 350 GeV due 
to the inclusion of QED effects. 
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Table 4: Evolution of the singlet part of the parton distribution in presence of QED 
effects from Qo = 2 GeV to Q = 100 GeV. 



The relative QED corrections, 



)2\ _ /withQED(Q ) fno QEd{Q 



r /'/^2\ J with C^EJJV^ J JnoytUV^ / { 

J no QEdI,V J 

in the parton distributions of the gluon, u- and d-quark are shown in Figures [T] as a 
function of x for Q = 100 GeV (l.h.s.) and as a function Q for x = 0.1 (r.h.s.). The 
QED effects are large and negative for large values of x, where the parton distributions 
are small. For x ^ 0.1, the QED corrections are < 0.3% and thus totally negligible. The 
r.h.s. of Figure n shows that the QED effects are negative and their size is increasing with 
increasing values of Q. As can be seen, they grow logarithmically with Q. 

In TablesHminiwe give the relative QED contributions of the evolved parton distribu- 
tions for several values of x and the evolution scales Q = 10 GeV, 100 GeV and 350 GeV. 
The contributions from QED corrections reach the order of 1% for large values of x. In 
Table ID we show results for the singlet part of the parton distributions for Q = 100 GeV 
and several values of x. All parton distribution functions decrease strongly with increasing 
X. As expected, the photon and distribution are rather small compared to the other 
distribution functions. However, we feel obliged to point out, that we made the ad- hoc 
assumption that the photon distribution 7 and the lepton distribution vanish at the 
input scale Qq. In a more realistic study the initial conditions have to be extracted from 
experimental data. 

4 Summary and conclusions 

We studied the inclusion of QED corrections into the evolution of parton distribu- 
tions. We showed that the usual set of parton distributions has to be enlarged by parton 
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distributions for the photon and for charged leptons. In addition, QED effects introduce 
additional mixing effects. To minimize the number of distributions, which mix under evo- 
lution, we introduced a new basis of parton distribution functions where only five partons 
mix. We reported on the implementation of the corrections into a numerical program 
and studied the QED effects in a toy model. The QED corrections can reach the order of 
1% and are large for large values of x. Hence, QED corrections to parton evolution are 
of the same order as expected contributions from NNLO corrections [15] and should be 
taken into account along with QCD NNLO corrections for precision observables, which 
significantly depend on large values of x. 
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